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Abstract 

We develop a method to compute the Casimir effect for arbitrary geometries. The 
method is based on the string-inspired worldline approach to quantum field theory 
and its numerical realization with Monte-Carlo techniques. Concentrating on Casimir 
forces between rigid bodies induced by a fluctuating scalar field, we test our method 
with the parallel-plate configuration. For the experimentally relevant sphere-plate 
configuration, we study curvature effects quantitatively and perform a comparison 
with the "proximity force approximation" , which is the standard approximation tech- 
nique. Sizable curvature effects are found for a distance-to-curvature-radius ratio of 
a/R> 0.02. Our method is embedded in renormalizable quantum field theory with 
a controlled treatment of the UV divergencies. As a technical by-product, we de- 
velop various efficient algorithms for generating closed-loop ensembles with Gaufiian 
distribution. 
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1 Introduction 



The Casimir effect has recently been under intense study, experimentally [2] as well as 
theoretically (for recent comprehensive reviews, see pj). In fact, we are currently witnessing 
a transition of the Casimir effect from a pure fundamental quantum effect, being interesting 
in its own right, via an experimentally challenging problem to a phenomenon becoming 
relevant to applied physics such as nanotechnology <A\. Moreover, the Casimir effect has 
been suggested as an experimentally powerful tool for investigating new physics beyond 
the standard model [5]. 

Considerable progress has been made in recent years as far as the Casimir effect of 
real (rather than idealized) conductors is concerned: the effects of finite conductivity, 
finite temperature, and surface roughness are theoretically well under control for the cur- 
rent experimental realizations. Even the dependence of the Casimir force on the isotopic 
composition of the interacting bodies has been studied recently By contrast, the depen- 
dence of the Casimir force on the geometry of the interacting bodies is neither completely 
understood nor quantitatively satisfactorily under control. Except for a small number 
of analytically solvable geometries, one has to rely on approximations among which the 
"proximity force approximation" |H] represents the most widely used method. Roughly 
speaking, the proximity force approximation maps the Casimir effect of an arbitrary geome- 
try onto Casimir's parallel-plate configuration, thereby neglecting curvature and tilt effects 
in an uncontrolled manner. In fact, the current limitations for a quantitative comparison 
of theory and experiment arise essentially from an estimated 1% error of the proximity 
force approximation. 

The basic obstacles against improving this situation are mainly technical in nature 
and partly fundamental. Standard strategies perform the Casimir calculations in two 
steps: first, the mode spectrum of quantum fluctuations in a given background geometry 
has to be identified; secondly, the Casimir energy is obtained by summing up (tracing 
over) the spectrum. The first step is obviously increasingly difficult the more complex a 
given geometry is; without a high degree of symmetry, even the use of standard numerical 
techniques is rather limited. The second step suffers from the same problems, but is 
moreover complicated by the fact that the mode sum is generally ultraviolet divergent. 
The divergencies have to be analyzed and, if possible, be removed by renormalization 
of physical parameters. Not only is the handling of these divergencies technically (and 
numerically) challenging, but the classification of divergencies is also still under intense 
debate 0HU1HI]. 

In this work, we propose a method that has the potential to solve these technical 
problems. Moreover, it is embedded in perturbative quantum field theory with its clear 
and unambiguous renormalization program. Our method is based on the "string-inspired" 
worldline formalism in which perturbative iV-point amplitudes are mapped onto quantum 
mechanical path integrals over closed worldlines ^2] (for a recent review, see ^H])- The 
technical advantages arise from the fact that the mode spectrum and its sum are not com- 
puted separately but all at once. These worldline integrals can conveniently be calculated 
with Monte-Carlo methods (worldline numerics) with an algorithm that is completely in- 
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dependent of the Casimir geometry; in particular, no background symmetry is required. 
Whereas the worldline integral is finite, the ultraviolet divergencies occur in a "propertime" 
integral, roughly corresponding to an integral over the size of the worldlines. The diver- 
gencies can be found at small propertimes (= small size = ultraviolet), where a mapping 
to Feynman-diagram language is possible and the standard rules of renormalization can be 
applied. 

In order to illustrate our method, we focus in this work on the calculation of Casimir 
forces between rigid bodies, induced by quantum fluctuations of a scalar field. The rigid 
bodies are modeled by background potentials V(x) (mainly of 5 function type), which allow 
us to approach the idealized limit of Dirichlet boundary conditions in a controlled way. As 
a benchmark test, we study the classic parallel-plate configuration in detail. Finally, we 
compute the Casimir forces between a plate and a cylinder as well as the experimentally 
highly relevant case of a plate and a sphere, both in the idealized Dirichlet limit. Here 
we find clear signals of curvature effects if the distance between the bodies is roughly a 
few percent of the cylinder /sphere radius or larger. This scale characterizes the limit of 
quantitative accuracy of the proximity force approximation. 

We developed the technique of worldline numerics in and it has successfully been 
applied to the computation of quantum energies or actions induced by scalar or fermion 
fluctuations in electromagnetic backgrounds [TH IT5| ITHj . As for any numerical method, 
possible finite-size or discretization errors have to be analyzed carefully. In this respect, 
the idealized Casimir problem turns out to be most challenging, because the background 
potentials with their £-like support affect the quantum fields on all scales. Therefore, we 
have to make sure that our worldline numerics operates sufficiently close to the "continuum 
limit" (propertime continuum in our case). We dedicate a whole section (Sect. Ej) to 
this question, also relevant for further applications of worldline numerics, and present a 
number of new and efficient algorithms for the generation of GauBian distributed closed- 
loop ensembles. 

Though the heart of our method is intrinsically numerical, we would like to emphasize 
that the worldline technique offers an intuitive approach to quantum phenomena. Par- 
ticularly for Casimir forces between rigid bodies, many features such as the sign of the 
interaction or curvature effects can easily be understood when thinking in terms of world- 
line ensembles (loop clouds). 

The paper is organized as follows: the next section provides a brief introduction into 
the worldline approach to the Casimir effect. Section El describes efficient methods for the 
generation of loop ensembles. The reader who is mainly interested in Casimir phenomenol- 
ogy may skip this section. Section 0] provides for an intuitive understanding of rigid-body 
Casimir forces in the light of the worldline language. Our numerical findings for the rigid- 
body Casimir force for several geometries (plate-plate, plate-sphere, plate-cylinder) are 
presented in section |S] 
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2 Worldline techniques for Casimir configurations 



2.1 Framework 

Let us discuss the formalism for the simplest case of a real scalar field <fi coupled to a 
background potential V(x) by which we describe the Casimir configuration. The field 
theoretic Lagrangian is 

C = \d^d^ + im 2 2 + hf{x) 4> 2 . (1) 

The potential V(x) can be considered as a spacetime dependent mass squared, implying 
that it has mass dimension 2. In the absence of any further fields and couplings, the 
complete unrenormalized quantum effective action for V is 

2 —a 1 + m z 



1 r ? / a 



2 ./i/A 2 ^ 



/ T U-T(-'9 2 W+V( a; ))| \ _ 1 r -m 2 T 

{ 1 1 ; (47rT)^ 2 



(3) 



Here we work in Z) = d + 1 Euclidean spacetime dimensions, i.e., d space dimensions. 
In Eq. (j3J), we have introduced the propertime representation of the Trln with UV cutoff 
A at the lower bound of the T integral. 1 Interpreting the matrix element as a quantum 
mechanical transition amplitude in propertime T, we can introduce the Feynman path 
integral, or worldline, representation, 

d D x (x\e- T (- d2+v ^\x) = J d D x CM M J Vxe -f T dT* 2 /*-£dTV( XCM+x (T))_ (4) 

x(0)=x(T) 

The Tr operation of Eq. (j2J), which has led to a transition amplitude at coincident points 
in Eq. 0, induces a path integral over closed worldlines, x(0) = x{T). In Eq. (|2J), we have 
shifted all worldline loops under the spacetime integral to have a common center of mass 
xcm, implying j Q drx^r) = 0. The normalization M is determined from the limit of zero 
potential, 

x(0)=x(T) 

such that the path integral can be interpreted as an expectation value with respect to an 
ensemble of worldlines with Gaufiian velocity distribution, 

\f f J) Tr -fndTX 2 /4-f^dTV(x CM +x(T))_ 1 / - J T cLt V (icm+^(t)) \ (a\ 

7 " (4vrT)^/ 2 \ h' [ } 

x(0)=x(T) 



^^Other regularization techniques are possible as well, e.g., dimensional regularization, (dT/T) — > 
^ 2€ (dT/r 1_e ); the propertime cutoff is used only for the sake of definiteness. For a pedagogical review of 
various regularization techniques in the Casimir context, see |17j . 
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The last step of our construction that is crucial for numerical efficiency consists of intro- 
ducing unit loops y(t) which are dimensionless closed worldlines parameterized by a unit 
propertime t G [0, 1], 



x^Tt) 



drx (t) 



dty 2 (t), 



(7) 



where the dot always denotes differentiation with respect to the argument. 

Inserting the path integral representation (jlj) into the effective action Q and using 
the unit loops y(t), we end up with the desired formula which is suitable for a numerical 
realization, 



r[v] = - 



i l 



dT 



-m 2 T 



2 (4tt) d / 2 J 1/A 2 T^+v/ 2 



d D x 



W v [y{t);x} 



Here and in the following we have dropped the subscript "CM" of the center-of-mass 
coordinate x M and introduced the "Wilson loop" 



W v [y(t);x] = exp 



— T / dtV(x + VTy{t)) 
Jo 



(9) 



and 



W v [y(t);x] 



J Vy W v [y{t);x]e-ti dt y 2 / A 
y(o)=y(i) 



(10) 



/ Vy e~ ft dt y 2/A 

y(0)=y(l) 

denotes the expectation value of an operator with respect to the path integral over unit 
loops y(t). This construction of Eq. (jSJ) is exact and completely analogous to the one 
proposed in jHj for electromagnetic backgrounds; further details can be found therein. 

For time-independent Casimir configurations, we can carry out the time integration 
trivially, J dxo = L xo , where L xo denotes the "volume" in time direction, and define the 
(unrenormalized) Casimir energy as 



£ — r/ L xo . 



2.2 Renormalization 

The analysis of divergencies in Casimir calculations is by no means trivial, as the ongoing 
debate in the literature demonstrates [TU1 ITT] . The reason is that divergencies in these 
problems can have different sources with different physical meaning. On the one hand, 
there are the standard field theoretic UV divergencies that can be mapped onto divergencies 
in a finite number of Feynman diagrams at a given loop order; only these divergencies can 
be removed by field theoretic renormalization, which is the subject of the present section. 

On the other hand, divergencies can arise from the modeling of the Casimir boundary 
conditions. In particular, idealized conditions such as perfectly conducting surfaces affect 
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quantum fluctuations of arbitrarily high frequency; therefore, an infinite amount of energy 
may be required to constrain a fluctuating field on all scales. These divergencies are real and 
imply that idealized conditions can be ill-defined in a strict sense. The physically important 
question is whether these divergencies affect the physical observable under consideration 
(such as Casimir forces) or not. If not, the idealized boundary conditions represent a 
simplifying and valid assumption, and the removal of these divergencies can be justified. 
But if the observable is affected, the idealized conditions have to be dropped, signaling the 
strong dependence of the result on the physical details of the boundary conditions (e.g., 
material properties). 

Even though the worldline is an appropriate tool for analyzing both types of divergen- 
cies, we concentrate on the first type in this paper, leaving a discussion of the second for 
future work. 

In order to isolate the field theoretic UV divergencies, we can expand the proper- 
time integrand for small propertimes (high momentum scales). Since this is equivalent 
to a local gradient expansion in terms of the potential V(x) (heat-kernel expansion), 
each term ~ V(x) n corresponds to a scalar one-loop Feynman diagram with n external 
legs coupling to the potential V(x) and its derivatives, and with the momentum inte- 
gration already performed (thanks to the worldline method). Using J dty^it) = and 
Jo dt (y^(t)y v {t)) y = (1/6)6^, we find up to order T 2 , 



which should be read together with the propertime factor \/T l+D ^ 2 in Eq. (jHJ). The term ~ 
V(x) corresponds to the tadpole graph. In the conventional "no-tadpole" renormalization 
scheme, the renormalization counter term ~ V(x) is chosen such that it cancels the tadpole 
contribution completely. Of course, any other renormalization scheme can be used as well. 
The corresponding counter term can be fixed unambiguously by an analysis of the tadpole 
Feynman diagram in the regularization at hand. In D < 4 spacetime dimensions, there is 
no further counter term, since V(x) has mass dimension 2. The remaining terms of 0(T 2 ) 
are UV finite in the limit T — > 1/A 2 — > 0. 

In 4 < D < 6, we need further subtractions. Here, it is useful to note that the last term 
on the first line of Eq. (J 12)) vanishes anyway, provided that the potential is localized or 
drops off sufficiently fast at infinity. This is, of course, always the case for physical Casimir 
configurations. 2 Renormalization provides us with a further counter term ~ J V 2 subject 
to a physically chosen renormalization condition such that the divergence arising from the 
last T 2 term is canceled. With this renormalization condition, the physical value of the 

2 Strictly speaking, infinitely extended surfaces such as idealized infinitely large plates do not belong to 
this class, but we can always think of large but finite surfaces and then take the infinite-surface limit after 
the infinite- volume limit. 




(12) 
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renormalized operator ~ V 2 is fixed. 3 For even higher dimensions, similar subtractions are 
required that involve higher-order terms not displayed in Eq. (|12j). 

As far as controlling divergencies by renormalization is concerned, this is all there is 
and no further ad hoc subtractions are permitted. However, having removed these UV 
divergencies with the appropriate counter terms does not guarantee that the resulting 
Casimir energy is finite. Further divergencies may arise from the form of the potential as 
is the case for the idealized Casimir energies mentioned above. 

In the present work, we take up a more practical position and are merely interested in 
the Casimir forces between disconnected rigid bodies which are represented by the potential 

V(x) = Vi(x) + V2(x) + We assume the rigid bodies as given, disregarding the problem 

of whether the Casimir energy of every single body is well defined by itself. For this, it 
suffices to study the interaction Casimir energy defined as the Casimir energy of the whole 
system minus the separate energies of the single components, 

E := £ v =Vi+v 2 +... — £v! — £v 2 — ■ ■ ■ ■ (13) 

Note that the subtractions do not contribute to the Casimir force which is obtained by 
differentiating the interaction energy with respect to parameters that characterize the 
separation and orientation of the bodies. By this differentiation, the subtractions drop 
out. Furthermore, these terms remove the field theoretic UV divergencies of Eq. (J 12)) : this 
is obvious for the terms linear in V(x); for the quadratic one, this follows from j V 2 = 
f (Vi + V 2 + . . . ) 2 = J (V 2 + V% + ■ ■ ■ ) • The last equation holds because of the local support 
of the disconnected bodies. By the same argument, the subtractions remove every term 
of a local expansion of 8v 1 +v 2 +... to any finite order. In this way, any divergence induced 
locally by the potentials is canceled. But, of course, the Casimir force is not removed - it 
is inherently nonlocal. 

The interaction energy in Eq. (fTSjl is also numerically favorable, since the subtractions 
can be carried out already on the level of the propertime integrands, avoiding manipulations 
with large numbers. 

We would like to stress that the definition of the interaction energy in Eq. (|13|) should 
not be confused with renormalization. It is a procedure for extracting exact information 
about the Casimir force between rigid bodies, circumventing the tedious question as to 
whether Casimir energy densities are locally well defined. This procedure also removes 
the field theoretic UV divergencies. In this case, renormalization conditions which fix the 
counter terms do not have to be specified. These local counter terms cannot exert an 
influence on the Casimir force for disconnected rigid bodies anyway, because the latter 
is a nonlocal phenomenon. Expressed in physical terms of the QED Casimir effect: the 
renormalized strength of the coupling between the electromagnetic field and the electrons 
in the metal is, of course, important for a computation of the local energy density near 

3 Since we used a gradient expansion, the renormalized operator is fixed in the small- momentum limit; if 
the renormalization condition operates at finite momentum, e.g., using the polarization operator, possible 
finite renormalization shifts can be obtained from an analysis of the corresponding Feynman diagram. 
However, in the present case of static Casimir problems, it is natural to impose a renormalization condition 
in the small-momentum limit anyway. 
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a plate, but the Casimir force between two plates is independent of the electromagnetic 
coupling constant. 

We would like to point out that the concept of the interaction energy is meaningless 
for the computation of Casimir stresses of single bodies, e.g., a sphere. Here, the renor- 
malization procedure has to be carried out as described above, and the result may depend 
on the renormalization conditions and strongly on the details of the potential. 



3 Worldline numerics 

In this section, we discuss possible numerical realizations of the worldline integral Eqs. (jHJ)- 
(I10J) (the more phenomenologically interested reader may proceed directly to Sect.EJ). 

As proposed in [Tl], we estimate the analytical integral over infinitely many closed 
worldlines by an ensemble average over finitely many closed loops obeying a Gaufiian 
velocity distribution P[{y(t)}], 

P[{y(t)}]=sU dtyitf) expf-~^W), with y(0) = y(l), (14) 

where the S constraint ensures that the loops are centered upon a common center of mass 
(here and in the following, we drop all normalizations of the distributions, because they 
are irrelevant when taking expectation values). Here, we have chosen to work with rescaled 
unit loopsy(t) as introduced in Eqs. (|S|)- (jl(J|) . Numerical arithmetics requires discretization; 
however, we generally do not discretize spacetime on a lattice, but only the loop propertime 
parameter t: 

{y(t)} - {y k }eR D , k = i,2,...,N, (15) 

where N denotes the number of points per loop (ppl). Whereas GauBian distributed 
numbers can easily be generated, the numerical difficulty is to impose the 5 constraint, 
Vx + V2 + • • • + Vn — 0, and the requirement of closeness. In the following, we discuss four 
possible algorithms, and recommend the last two of them based on Fourier decomposition 
("f loops") or explicit diagonalization ("v loops"). 



3.1 Heat-bath algorithm 

A standard approach for the generation of field (or path) distributions that obey a certain 
action is the heat-bath algorithm, which has been employed for worldline numerics in 
[HJ [T3J HH]. Discretizing the derivative in the exponent of Eq. dHJ), e.g., by y — > N{yk — 
yk-i), each point on a loop can be regarded as exposed to a "heat bath" of all neighboring 
points. The discretized probability distribution then reads 

N N 

P[{Vk}} = s(y! + ... + y N ) expl--^^-^) 2 } . (16) 

fc=i 
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where yo = Vn- The heat-bath procedure now consists in the following steps: (i) choose a 
site i G [l,iV], consider all variables y^, k i as constant, and generate the y^ according 
to its probability; (ii) visit all variables of the loops (e.g., in a serial fashion or using the 
checkerboard algorithm). Thereby, the closeness requirement is easily realized with, e.g., 
yN being in the heat bath of i/n~i and y\, etc. The center-of-mass constraint can be 
accommodated by shifting the whole loop correspondingly after one thermalization sweep 
(update of all points per loop). 

Whereas this procedure has been sufficient for the applications discussed in [TJJ [Tj)J [TE] , 
it turns out that this algorithm suffers in practice from a thermalization problem for large 
values N. To demonstrate this, let us define the extension e of the loop ensemble by the 
loop mean square 

e 2 = ^ J d yi ...dy N y 2 k P[{y k }] , M = J d Vl . . . dy N P[{y k }] . (17) 
This quantity can be calculated analytically, straightforwardly yielding 

e = fii 1 -^)- (18) 

In order to generate a "thermalized" loop, one starts with a random ensemble {y k } and 
performs n t heat-bath sweeps. For each loop, we calculate its extension e. After averaging 
over 1000 loops, we compare the estimator of e as function of n t with the analytic result 
(JTHJ) corresponding to the limit n t — > oo. 

The result is shown in Fig. ^ One clearly observes that the thermalization of loop 
ensembles is expensive for iV > 500. In fact, roughly n t = 45000 is needed for an acceptable 
loop ensemble consisting of iV = 1000 points. Since a computation of Casimir energies 
requires loop ensembles of N > 1000, the heat-bath algorithm is too inefficient and cannot 
be recommended. 



3.2 Random Walk 

In order to circumvent the thermalization problem, one may exploit the connection between 
loops with GauBian velocity distribution and random walks [THIHH]. This has been adapted 
to worldline numerics with latticized spacetime in [21]; here, however, we keep spacetime 
continuous. For this purpose, let us give up the concept of unit loops for a moment, and 
reinstate the naturally emerging coordinate space loops x(t), 

x(t) = 4=v(r/T), x(0)=x(T). (19) 

Probability theory tells us that random walks automatically implement the GauBian 
velocity distribution 

N-l , x . 

ri exp \~AK^ {xi+i ~ xi)2 \- (2o) 

1=1 ^ ' 
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Figure 1: The average extension e (multiplied by y N for better visualization) of the loops 
as function of the number of thermalizations n t . 



The crucial point is to establish the relation between a loop that a random walker with 
step length s would generate for us and a thermalized loop at a given propertime T. This 
relation results from a coarse-graining procedure, which we present here briefly. Given that 
the random walker starts at the point Xi, the probability density for reaching the point Xf 
after n steps is given by 

f n_1 1 

p(xf | Xi, n,s) = J d D x 2 • • • d D x n _i JJ q^ s d-i S ^\ Xk+1 ~ Xk ' 

with Q(D) being the solid angle in D dimensions, X\ = Xi and x n = Xf. For n»l, but 
ns 2 fixed, the central-limit theorem can be applied |19j : 



ns 



fixed. 



Comparing (|2*T|l with (|20|) . one identifies 



At 



ns 
2D 



(21) 



(22) 



The dimension of the propertime as well as its relation to the loop length L appears here 
in an obvious way, 

■ V "" s2 LS (23) 



T 



2D 2D' 

where iV\„ now is the total number of walker steps. Is is important to point out that 
the propertime can be tuned in two ways: we can adjust either the walker step s or the 
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number of points N w . The corresponding two methods to generate a loop ensemble at 
given propertime T work as follows: 

Method 1 : s is fixed. 

(1) choose the walker step s; 

(2) read off from Eq. (}2*3*j) the number of points N w corresponding to T; 

(3a) generate N w points by letting a random walker go N w steps, and accept the config- 
uration if the last step leads him into a small sphere (radius e) centered upon the 
starting point; 

(3b) close the loop 'by hand' by shifting the last point to the starting point; 

(4) shift the center of mass to zero; 

(5) repeat steps (3) and (4) times for an ensemble of loops. 

We point out that the value of s must be much smaller than the characteristic length 
scale provided by the background potential. A second constraint on s arises from the 
applicability of the central limit theorem, i.e., n ^> 1 in ([21)1 . A third systematic numerical 
uncertainty follows from the shift in step (3b). Unfortunately, small values for e result in 
low acceptance rate for loops, and, therefore, increase the numerical effort to generate the 
loop ensemble. A good compromise is to set the radius e to some percentage of the step 
length s. 

For illustration only, we leave the Casimir effect for a second and consider the average 
Wilson Loop (Wy) (see Eq. (jIDj) ) for the case of a constant magnetic background field 
B = Be z at T = 1 and D = 2, 

V(x) = A k (x) x k , A = B/2 (y, -x, 0) . 

For T = 1 the walker step length is given by s = -7^=. Figure El shows our numerical result 
as a function of N w in comparison with the exact value. Circles with error bars correspond 
to loop ensembles generated with e = 0.05 s. The limit ()21|) seems to be attained for 
50 < N w < 100 (s < 0.3). For a further improvement of the numerical accuracy, large 
values of N w and a decrease of e at the same time are required. Finally, we point out that 
the deviation from the exact result in the case of the heat-bath-generated loop ensemble 
(blue square) is probably due to thermalization effects. 

Note that we have to generate a loop ensemble for each value of T (~ N w s 2 ), which 
makes this procedure far more memory consuming than the heat-bath approach. If we 
decide to generate the loop ensembles once and for all and save them to disk, we have to 
handle huge amounts of data. On the other hand, if we create our loops 'on demand' (while 
performing the T or x integrations), we are confronted with a serious waste of computing 
time. 

Method 2 : N w is fixed. 
(1) choose the number of points iV^,; 
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Figure 2: Average Wilson Loop (Wy) (cf. Eq. (I10JI ) for the case of a constant magnetic 
background field B for B = 1, T = 1 and D = 2 as a function of the number of points per 
loop. 



(2) set the walker step to s — 1; 

(3) proceed with steps (3), (4) and (5) of the first method. 

The loop ensemble is here generated only once and then rescaled to adjust the step length 
to the value s corresponding to T in (|23|) . This method therefore works as in the case of 
the rescalable thermalized unit loops, with the difference that the propertime rescaling is 
realized via the rescaling of the step length. This tuning at the level of s provides for a 
better control of the microscopic features of the loops. The second procedure is thus a good 
candidate to replace the thermalized loops since it combines the absence of thermalization 
and the rescaling of an all-at-once generated ensemble. 

It should however be emphasized that most of the computer time is spent on generating 
redundant open loops. This is due to the fact that, for a given e, the fraction of loops 
which close after N w steps decreases like Nw D ^ 2 - 



3.3 Fourier decomposition: "f loops" 



We are now looking for alternative methods that could combine some advantages of the 
two previous approaches and bypass the problems rendering them impractical. A highly 
efficient procedure arises from a Fourier representation of our unit loops 



N 



y(t) = E 



u=0 



a,, cos 



( 2tt, 



vt ) + b u sin ( 2n v t 



a 







(24) 



where N is the number of Fourier modes included (which agrees with the number of points 
specifying each loop, see below). The choice a = guarantees that the loop center of 
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mass is located at the origin. Inserting Eq. into Eq. (JHj), the probability distribution 
for the coefficients is given by 

2 N 2 N 

P[a,b] = exp{-^-$>2 - y£^}- (25) 

u=l u=0 

We can then take advantage of the fact that the Fourier components {a, b} are not cor- 
related, in order to generate our loops in momentum space. The reconstruction of the unit 
loop y(t) in Eq. (J24|) is most efficiently performed by using the fast Fourier transformation 
(FFT). For these purposes, we define complex coefficients c u := a u + ib u , and obtain 

y (t) = ^^ Cu exp|-z27rtz/| . (26) 

v 

The FFT procedure generates a series of points yi, i = ... iV — 1 which discretize the 
continuous curve y(t) and thereby constitute the unit loops. 



3.4 Explicit diagonalization: "v loops" 

Finally, we propose an algorithm that is based on a linear variable transformation {y^} — > 
{vk}, such that the discretized distribution (jTUJl becomes purely GauBian. These new 
variables are velocity- like and diagonalize the quadratic form in the exponent. 

Because of the 6 function in Eq. (|16jh only N — 1 coordinates per loop are independent. 

N 

poo 



Defining J Vy = FJ dyi, we may perform, e.g., the yjy integration using the 5 function, 

00 i=i 

r r N ~ x 

I VyP[{y k }] ... = I \\_dyi e [-f(E£2 1 (^-K-i) 2 +( 2 s/i+y2+-+s/iv-i) 2 +(!/i+y2+-+2 2/J v„ 1 ) 2 )] 
J i=i 

N-l 

II . . . , (27) 



1=1 

where the dots represent an arbitrary y-dependent operator, and we introduced the abbre- 
viation Y for the quadratic form. In order to turn the exponential into a product of simple 
GauBians, we define N — 1 new velocity-like variables, 

3 3 
Vi := + 2/2 + 2/3 H h Vn-2 + ^Vn-i, 

Vi := yi-yi-!, i = 2,3,...,N-l. (28) 
For notational simplicity, it is useful to also introduce the auxiliary variable, 

Vi,j = Vi + Vi-x -\ \- v j+1 = yi - yj, for % > j = 1, 2, . . . , N - 1, (29) 
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such that the exponent Y can be written as 

N-l 



Y 



i=2 



1 



J2 V i + ( V l - ) + ( V l + y V N-l,l 



N-l 



2v 



2 , 1 2 i V"^ 2 

i + ^jv-i,i + • 



(30) 



i=2 



We observe that the variable vi now appears quadratically in the exponent as desired. The 
same has still to be achieved for v 2 ■ ■ ■ ^jv-i- For this, we note that ujv-i,i — v N-i + v n-2,i 
by definition (|29|). Defining 



1 



(31) 



we indeed obtain for the exponent Y 

Y = 2v 2 1 +v 2 N _ 1 + ^(v N _ 1 +v N ^2,i) 2 + J2 
3 1 N ~ 2 



JV-2 



i=2 



(32) 



i=2 



where vjj_ 1 also appears quadratically. We can continue this construction by defining 

1 



vn- 



v N-i + — —x v N-i-i,h i = l,...,N -2 , 



i + 2 



which turns the exponent Y into a purely GauBian form: 

Y = 2v\ + \v\_ x + \v 2 N _ 2 + ■■■ + ^ift-i 
2 3 i + 1 



H h 



iV 



iV- 1 



(33) 



(34) 



The last step of this construction consists in noting that we can substitute the integration 
variables according to 



N-l 



N-l 



N-l 



Y[ dyi = J \\ dvidvi = J Y\ dvi = Vi 



(35) 



i=i 



i=2 



1=1 



with nonzero but constant Jacobians J, J, the value of which is unimportant for the 
calculation of expectation values. This allows us to write the path integral Eq. (j2ZJ) as 



J VyP[{y k }]--- = J Jvv exp 



N 

T 



JV-2 



i=l 



= J / VvP[{v k }} 



(36) 

where P[{vfc}] can now be generated straightforwardly with the Box-Miiller method |2()j . 

For the construction of unit loops ("v loops"), the above steps have to be performed 
backwards. The recipe is the following: 
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(1) generate N — 1 numbers Wi, i — 1, . . . , N — 1 via the Box-Miiller method such that 
they are distributed according to exp(— wf); 

(2) compute the i — 1, . . . , TV — 1, by normalizing the wf 

2 



''J = \IJ^ W ^ 



(3) compute the Vi, i — 2, . . . , N — 1, using 

1 4-1 

''' = ! ~ N + 2 - ■ Vi ~ 1 ' 1 ' where ''' :-i X! : ( 38 ) 

* J=2 

(4) construct the unit loops according to 



" 5('-i) 



Vi = yi-i + Vi, i = 2, . . . ,N - 1, 

JV-l 

= - yi ; (39) 
1=1 

(5) repeat this procedure «l times for ul unit loops. 

The formulas in step (4) can be checked straightforwardly by inserting the definitions 
of the ViS and v±. 

This v-loop algorithm allows us to generate unit loops efficiently without thermalization, 
i.e., no redundant thermalization sweeps have to be performed, and works for an arbitrary 
number of points per loop N. 



3.5 Benchmark test 

We test the quality of our loops with the aid of the Casimir energy for the parallel-plate 
configuration in the Dirichlet limit, the physics of which is described in the next section. 

As far as numerics is concerned, there are basically two parameters that control the 
quality of our loop ensemble: the number of points per loop (ppl) N, and the number of 
loops n L . The larger these numbers, the more accurate is our numerical estimate at the 
expense of CPU time and size. Whereas increasing the number of loops reduces the 
statistical error of the Monte-Carlo procedure, increasing the number of ppl iV reduces the 
systematic error of loop discretization. 



15 



-1.7 



-1.9 



pa 



-2 — 



-2.1 



-2.2 



I 



i 



~r 



□ 100 000 f loops 
A 1500 v loops 
O 100 000 v loops 
exact result 



_L 



1000 10000 le+05 

N (points per loop) 



-2 <g 



le+06 



Figure 3: Numerical estimate of the interaction Casimir energy of the parallel-plate config- 
uration for various loop ensembles as a function of the number of points per loop N. The 
error bars correspond to the Monte-Carlo statistical error; deviations from the exact result 
on top of the statistical error measure the systematic error due to loop discretization. 

In order to estimate this systematic error, we have to study the approach towards the 
continuum limit. The idea is to choose N large enough for a given nL, such that the 
systematic error is smaller than the statistical one. 

In Fig. |3J we plot the numerical estimates for the parallel-plate Casimir energy as a 
function of the number of ppl N and compare it with the classic result. The error bars 
represent the statistical error of the Monte-Carlo procedure. The deviation of the numerical 
estimates from the exact result on top of the error bars serves as a measure of the systematic 
error. As is visible therein, a rather small number of several thousand ppl, iV > (9(1000), 
is sufficient to get a numerical estimate with < 5% error using = 1500 loops. For a 
high-precision estimate with an error < 0.5%, larger loop ensembles with > 100 000 are 
required. For N ~ 50 OOOppl, systematic and statistical errors are of the same order, and 
for iV > 100 OOOppl, the systematic error is no longer relevant for v loops. For f loops, 
however, we observe a systematic 1% error in the high-precision data of unclear origin. 

Nevertheless, the important conclusion of this test is that worldline numerics has proved 
its ability to describe quantum fluctuations with Dirichlet boundary conditions quantita- 
tively. 
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4 Casimir forces between rigid bodies 



Casimir forces can be analytically computed for only a small number of rigid-body geome- 
tries among which there is Casimir's classic result for the parallel-plate configuration; for 
perfectly conducting plates at a distance a, the interaction energy per unit area is [I] 



for a fluctuating real scalar field; for a complex scalar as well as for electromagnetic fluc- 
tuations, the factor 1/2 has to be dropped. The famous Casimir force is obtained by 
differentiating Eq. (jjOj) by a. 

4.1 Proximity force approximation 

The standard approximation method for not analytically solvable Casimir problems is the 
proximity force approximation (PEA) IB]- The basic idea is to apply the parallel-plate 
result to infinitesimal bits of the generally curved surfaces and integrate them up, 



where Epp is the interaction energy per unit area of the parallel-plate case. S represents 
the integration domain and denotes either one of the surfaces of the interacting bodies or 
a suitably chosen mean surface [H]. At this point, the proximity force approximation is 
ambiguous, and we will simply insert both surfaces in order to determine the variance. 
In Eq. ([41)1 . da denotes the invariant surface measure, and z represents the separation 
between the two surfaces associated with the surface element da on S. Obviously, the 
proximity force approximation neglects any nonparallelity and any curvature - the latter 
because each surface element on SI is assumed to "see" only one surface element on S2 at 
separation z\ but curvature effects require information about a whole neighborhood of the 
element on S2. 

The proximity force approximation is expected to give reasonable results only if (i) the 
typical curvature radii of the surfaces elements is large compared to the element distance 
and (ii) the surface elements with strong nonparallelity are further separated than the more 
parallel ones. 4 

For configurations that do not meet the validity criteria of the proximity force approx- 
imation, a number of further approximations or improvements exist, such as an additive 
summation of interatomic pairwise interactions and the inclusion of screening effects of 
more distant layers by closer ones [3] . Though these methods have proved useful and even 
quantitatively precise for a number of examples, to our knowledge, a general, unambiguous 
and systematically improvable recipe without ad hoc assumptions is still missing. 

4 The second condition is not so well discussed in the literature; it is the reason why the proximity force 
approximation gives reasonable results for a convex spherical lens over a plate (convex as seen from the 
plate), but fails for a concave lens. 




(40) 




(41) 
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fa) (b) (c) 



Figure 4: Worldline loop contributions to Casimir energies between two surfaces (SI and 
S2): loop (a) does not contribute at all, it is an ordinary vacuum fluctuation. Loop (b) 
contributes to the local energy density near the upper plate, but does not contribute to the 
Casimir force. Only loop (c) contributes to the Casimir force, since it "sees" both surfaces. 
Here, the loop picks up nonlocal information about a whole neighborhood, whereas the 
proximity force approximation employs only information about local distances indicated 
by the dashed line. 

In Sect. we compare our results with the proximity force approximation in the sim- 
plest version as mentioned above, in order to gain insight into the effects of curvature. 

4.2 Casimir forces on the worldline 

As described in Sect. 121 we represent the rigid bodies by a potential V(x). The functional 
form of the potential leaves room enough for modeling many physical properties of real 
Casimir configurations. Let us confine ourselves to an idealized potential well which is 
represented by a 5 function in space (for "soft" boundary conditions, see, e.g., |2*2*j). 



where the geometry of the Casimir configuration is represented by S, denoting a d — 1 
dimensional surface. £ is generally disconnected (e.g., two disconnected plates, £ = Si+S 2 ) 
and can be degenerate, i.e., effectively lower dimensional (a point). The surface measure 
da is assumed to be reparametrization invariant, and x a denotes a vector pointing onto 
the surface. The coupling A has mass dimension 1 and is assumed to be positive. It can 
roughly be viewed as a plasma frequency of the boundary matter: for fluctuations with 
frequency uo ^> A, the Casimir boundaries become transparent. In the limit A — > oo, the 
potential imposes the Dirichlet boundary condition, implying that all modes of the field 
have to vanish on E. 




(42) 
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Inserting this potential into the worldline formula (JBJ), we encounter the integral 




da Z («) 

where {ij} is the set of all points where a given scaled unit loop y/Ty(t) centered upon x 
pierces the Casimir surface S at If a loop does not pierce the surface (for given T and 
x), Iy[y(t)] = for this loop. Of course, there are also loops that merely touch the Casimir 
surface but do not pierce it. For these loops, the inverse velocity l/\y(ti)\ diverges on the 
surface. But since this divergence occurs in the argument of an exponential function, these 
loops remove themselves from the ensemble average. 

As an example, let £ consist of two disconnected surfaces (bodies), such that V(x) = 
Vi(x) + V2(x). For a given propertime T, the Casimir energy density at point x receives 
contributions only from those loops which pierces one of the surfaces. The interaction 
energy density defined in Eq. (I1HJ) is even more restrictive: if a certain loop yo(t) does 
not pierce one of the surfaces, then (Wv 1+ v 2 [j/o] — 1) — (Wyjz/o] — 1) — (WV 2 [2/o] — 1) = 0. 
Therefore, only those loops which pierce both surfaces contribute to the interaction energy 
density, as illustrated in Fig. HJ 

If the loop yo(t) does pierce both plates, its contribution to the energy density is 

contrib. of y Q (t) = 1 - (e _T/v i + e~ TIv ^ - e~ TIv ^) £ (0, 1]. (44) 

From this general consideration, together with the global minus sign in Eq. (jHJ), we learn 
that Casimir forces between rigid bodies in our scalar model are always attractive. This 
statement holds, independent of the shape of the bodies and the details of the potential 
(as long as V(x) is non-negative). 

In the Dirichlet limit, A — > oo, the exponential functions in Eq. ()44|) vanish, and the 
contribution of a loop is = 1 if it pierces both surfaces and = otherwise. 

5 Numerical results 
5.1 Parallel Plates 

Let us first consider the classic example of a Casimir configuration consisting of parallel 
plates separated by a distance a and located at z — —a/2 and z = a/2 orthogonal to the 
z = Xd axis. For this, Eq. reduces to 

V(x) = V(z) = X[8(z + a/2) + 5(z - a/2)} =V X + V 2 . (45) 

In order to test the numerical worldline approach, we compare our numerical estimates 
with the analytically known result [2H1 of the interaction Casimir energy for arbitrary 
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X a 

Figure 5: Parallel plates: interaction Casimir energy per unit area for the parallel-plate 
configuration as a function of the coupling A (units are set by the plate separation a) . The 
numerical estimate reproduces the exact result for a wide range of couplings including the 
Dirichlet limit (cf. FigJ2J). 



coupling A and scalar mass m in units of the plate separation a. In Fig. we study a wide 
range of couplings and the approach to the Dirichlet limit, Aa ^> 1; here, the energy per 
unit area tends to 

n „ , 1 7T 4 1 1 1 

hm Ep P {\, a) = — — — ~ — — — x 2.16 . . . x — , 46 

Aa^oo FFK ' ' 2(4tt) 2 45 a 3 2(4vr) 2 a 3 ' v ; 

which is the classic Casimir result for a massless scalar field. 5 As is visible in Fig. EJ the 
agreement is satisfactory even for small ensembles with N = 20'OOOppl. 

Let us finally discuss the Casimir energy as function of the distance a of two parallel 
plates for finite mass m and finite A, in order to explore the strength of the worldline 
approach in various parameter ranges. The result is shown in Fig. |3 A finite value 
for A simulates a finite plasma frequency. Hence, for a C 1/A the plates become more 
transparent for those modes of the quantum field which fit between the plates. This weakens 
the increase of the interaction Casimir energy for decreasing plate separation, which turns 
from ~ 1/a 3 into a A 2 /a law 23J. For a ^> l/m, we observe that the Casimir energy 
decreases exponentially with a, as expected, since possible fluctuations are suppressed 
by the mass gap. In the intermediate distance regime, 1/A a < l/m, a reasonable 

5 Here and in the following, we have explicitly displayed the common propertime prefactors l/[2(47r) £) / 2 ] 
for convenience (see prefactor in Eq. |JS}). 
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Figure 6: Parallel plates: interaction Casimir energy per unit area for the parallel-plate 
configuration as a function of the distance a in units of the mass m for A = 100m. The 
exact result (solid line) j2Sj is well reproduced by the numerical estimate over many orders 
of magnitude. For intermediate parameter values, the classic Casimir result (idealized 
Dirichlet limit Eq. (|46|h dashed line) represents a reasonable approximation. 

approximation is given by the classic power law Epp ~ 1/a 3 , which is familiar from the 
ideal case A — > oo, m = 0. 

5.2 Sphere above plate 

The Casimir force between a sphere or a spherical lens above a plate is of utmost impor- 
tance, because a number of high-precision measurements have been performed with this 
experimental configuration. Let us confine ourselves to the massless case, m = 0, in the 
Dirichlet limit A — > oo; generalizations to other parameter ranges are straightforward, as 
in the parallel-plate case. 

In order to gain some intuition for curvature effects, let us consider a sphere of radius R 
the center of which resides over a plate at distance a = R as an example. The interaction 
Casimir energy density along the symmetry axis is shown in Fig. For comparison, the 
energy density of the case where the sphere is replaced by a plate is also shown. One 
observes that the energy density close to the sphere is well approximated by the energy 
density provided by the parallel-plates scenario. This is already at the heart of the nonlocal 
nature of the Casimir force and can easily be understood in the worldline approach. 
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Figure 7: Sphere above Plate: interaction Casimir energy density along the symmetry 
axis (x axis) for the sphere-plate configuration in comparison to the parallel-plate case. 
Close to the sphere, the worldline loops do not "see" the curvature; but at larger distances, 
curvature effects enter the energy density. For illustration, the sphere-plate geometry is 
also sketched (thin black lines). 

Recall that the dominant contribution to the interaction Casimir energy density arises 
from loops which intersect both surfaces. If the center of the loop is located close to 
the sphere, the loops which intersect both surfaces hardly experience the curvature of 
the sphere; this is because loops that are large enough to pierce the distant plate will 
also pierce the close-by sphere rather independent of its radius. By contrast, if the loop 
center is located close to the plate, the dominant (large) loops possess intersections with 
the sphere at many different points - not necessarily the closest point. In this case, the 
worldline loops "see" the curvature of the sphere that now enters the energy density. 

Let us now consider the complete interaction Casimir energy for the sphere-plate config- 
uration as a function of the sphere-plate distance a (we express all dimensionful quantities 
as a function of the sphere radius R). In Fig. |HJ we plot our numerical results in the range 
a/R ~ (9(0.001 . . . 10). Since the energy varies over a wide range of scales, already small 
loop ensembles with rather large errors suffice for a satisfactory estimate (the error bars of 
an ensemble of 1500 v loops with 4000 ppl cannot be resolved in Fig. [HJ). 

Let us compare our numerical estimate with the proximity force approximation (PFA): 
using the plate surface as the integration domain in Eq. (JHJ), S = S^te, we obtain the 
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Figure 8: Sphere above Plate: logarithmic plot of the interaction Casimir energy for the 
sphere-plate configuration. For small separations/large spheres, a/R < 0.02, the proximity 
force approximation (PFA) approximates the numerical estimate well; but for larger a/R, 
curvature effects are not properly taken into account. The PFA becomes ambiguous for 
larger a/R, owing to possible different choices of the integration domain S in Eq. (JHJ. 
A geometric mean (dotted-dashed line) of S — S p i a t e and S = Sphere shows reasonable 
agreement with the numerical result. 

solid line in Fig. |H] (PFA, plate-based), corresponding to a "no-curvature" approximation. 
As expected, the PFA approximation agrees with our numerical result for small distances 
(large sphere radius). Sizable deviations from the PFA approximation of the order of a few 
percent occur for a/R~ 0.02 and larger. Here, the curvature-neglecting approximations 
are clearly no longer valid. This can be read off from Fig. El where the resulting interaction 
energies are normalized to the numerical result. 

In the PFA, we have the freedom to choose alternatively the sphere surface as the 
integration domain, S = Sphere- Although still no curvature-related fluctuation effects 
enter this approximation, one may argue that information about the curvature is accounted 
for by the fact that the integration domain now is a curved manifold. Indeed, Fig. |H] shows 
that this "sphere-based" PFA approximation deviates from the plate-based PFA in the 
same direction as the numerical estimate, but overshoots the latter by far. It is interesting 
to observe that the geometric mean, contrary to the arithmetic mean, of the two different 
PFA approximations lies rather close to the numerical estimate; we will comment on this 
in more detail in the next section. 
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Figure 9: Sphere above Plate: interaction Casimir energies normalized to the numerical 
result (further conventions as in Fig.|SJ). For a/R> 0.02, the fluctuation-induced curvature 
effects occur at the percent level. 

5.3 Cylinder above plate 

In order to study the relation between PFA approximations and the full numerical estimate 
a bit further, let us consider a second example of a cylinder above a plate. Apart from the 
difference in the third dimension, all parameters and conventions are as before. 

Again, we observe in Fig. El that the numerical estimate is well approximated by 
the PFA for a/R < 0.02, but curvature effects become important for larger distance-to- 
curvature-radius ratios. As in the sphere-plate case, the plate-based PFA neglects, but the 
cylinder-based PFA over-estimates, the curvature effects for a/R of order one. 

Our results seem to suggest that the various possible choices for the integration domain 
in the proximity force approximation may give upper and lower bounds for the correct 
answer. Indeed, the geometric mean between the two possible choices for the sphere- 
plate configuration is rather close to the numerical estimate (dotted-dashed line in Figs. El 
and ITn|) . Similar positive results for the geometric mean have been found for the two- 
concentric-cylinder configuration [21] using semiclassical approximations and for a 
"chaotic" geometry 

However, we believe that this "agreement" beyond the strict validity limit of the PFA 
is accidental. First, detailed inspection reveals that the geometric mean and the numeri- 
cal estimate are not fully compatible within error bars; this is particularly visible in the 
cylinder-plate case in Fig. Secondly, there are no fundamental arguments favoring the 
geometric mean; by contrast, the arithmetic mean (as well as the quadratic mean) are not 
good approximations. Thirdly, for even larger separations, a/R —>■ oo, it is known that the 
interaction Casimir energy in the sphere-plate case behaves as ~ i? 3 /a 4 |2b| . whereas even 
the sphere-based PFA decreases only with ~ R 2 /a 3 . From the viewpoint of the worldline, 
it is obvious anyway that true fluctuation-induced curvature effects cannot be taken into 
account by PFA-like arguments. Nevertheless, the geometric-mean prescription may yield 
a reasonable first guess for Casimir forces in a parameter range beyond the formal validity 
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Figure 10: Cylinder above Plate: logarithmic plot of the interaction Casimir energy for 
the sphere-plate configuration (cf. Fig. [SJ). 



bounds of the PFA where the expansion parameter is maximally of order one. 



6 Conclusions 

We have proposed and developed a new method to compute Casimir energies for arbitrary 
geometries from first principles in a systematic manner. The approach is based on pertur- 
bative quantum field theory in the string-inspired worldline formulation which maps field 
theoretic problems onto one-dimensional quantum mechanical path integrals with an evo- 
lution in a "5th coordinate" , the propertime. These path integrals can easily be performed 
with numerical Monte-Carlo techniques. 

Beyond any technical and numerical advantages, we first would like to stress that the 
worldline formulation offers an intuitive approach to the phenomena induced by quantum 
fluctuations. The geometric dependence of Casimir forces between rigid bodies, curvature 
effects and nonlocalities can already be guessed when thinking in terms of worldline loop 
clouds. 

As to technical advantages, the (usually complicated) analysis of the fluctuation spec- 
trum and the mode summation are performed at one fell swoop in the worldline approach. 
Above all, our algorithm is completely independent of the details of the Casimir geometry 
and no underlying symmetry is required. The algorithm is scalable: if higher precision is 
required, only the parameters of the loop ensemble (points per loop and number of loops) 
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have to be adjusted 6 . 

In this work, we have focused on Casimir forces between rigid bodies for which a 
computation of the interaction energy suffices; the latter is free of subtle problems with 
renormalization. Nevertheless, the worldline approach is in principle capable of isolating 
and classifying divergencies of general Casimir energy calculations, and the unambiguous 
program of quantum field theoretic renormalization can be performed. 

Confining ourselves to a fluctuating real scalar field, we tested our method using the 
parallel-plate configuration. New results have been obtained for the experimentally im- 
portant sphere-plate configuration: here we studied the (usually neglected) nonlocal cur- 
vature effects which become sizable for a distance-to-curvature-radius ratio of a/R > 0.02. 
Even though the proximity force approximation (PFA) as standard approximation method 
cannot correctly account for fluctuation-induced curvature effects, we found (accidental) 
agreement between our numerical estimate and the PFA with a "geometric-mean prescrip- 
tion": the latter implies a geometric mean over the possible choices of surface integration 
in Eq. (|41j). This geometric mean PFA might provide for a first guess of the Casimir force 
for a/R of order one, but has to be treated with strong reservations. 

In this work, we have accepted a number of simplifications, in order to illustrate our 
method. Many generalizations to more realistic systems are straightforward, as discussed 
in the remainder of this section: 

1) We modeled the Casimir bodies by 6 potentials, mostly taking the Dirichlet limit. In 
fact, this was not a real simplification, but numerically even more demanding. Modeling 
the bodies by finite and smooth potential wells requires worldline ensembles with a much 
smaller number of points per loop. The 5 potentials represent the "worst case" for our 
algorithm, which has nevertheless proved to be applicable. 

2) In experimental realizations, effects of finite temperature and surface roughness have to 
be taken into account. Both can be implemented in our formalism from first principles. 
Including finite temperature with the Matsubara formalism leads to a worldline integral 
with periodic boundary conditions of the worldline loops in Euclidean time direction [2*7| ITo] 
which can easily be performed for Casimir configurations. The surface roughness can 
be accounted for by adding a characteristic random "noise" to the local support of the 
potential. In both cases, the observables can directly be computed by our formalism 
without any kind of perturbative expansion. 

3) For obtaining the Casimir force, our results for the interaction energy have to be differen- 
tiated with respect to the separation parameter. Since numerical differentiation generally 
leads to accuracy reduction, it is alternatively possible to perform the differentiation first 
analytically; this yields a slightly more complicated worldline integrand which can never- 
theless be easily evaluated without loss of precision. By a similar reasoning, we can also 
obtain the (expectation value of the) energy-momentum tensor, which is frequently at the 
center of interest in Casimir calculations. For this, we can exploit the fact that the energy- 

6 The numerical computations for this work have been performed on ordinary desktop PC's. Improve- 
ment in precision can be obtained at comparatively low cost, since the computer resources required increase 
only linearly with our loop parameters. 
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momentum tensor can be obtained from the effective action by differentiating Eq. (|SJ) with 
respect to the metric analytically; the resulting worldline integrand can then be put into 
the standard path integral machinery. 

4) Radiative corrections to the Casimir effect can also be included in our method, employing 
the higher- loop techniques of the worldline approach We expect these computations 
to be numerically more demanding, since more integrations are necessary, but the general 
approach remains the same. 

5) The implementation of finite conductivity corrections is less straightforward, since this 
generally requires a formulation for real electromagnetic fluctuations (an extension to com- 
plex scalars is not sufficient). For this, the starting point can be a field theoretic Lagrangian 
defining a model for the interaction of the electromagnetic field with the bodies as sug- 
gested, e.g., in [2H]- Although these Lagrangians are generally not renormalizable, one may 
expect that the dispersive properties of the bodies provide for a physical ultraviolet cutoff 
(although this has to be studied with great care 29 J. 
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